Charged particle (pseudo-) rapidity distributions in p+p/p+p and Pb+Pb/Au+Au 
collisions from SPS to LHC energies from UrQMD 
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We present results for final state charged particle (pseudo-)rapidity distributions in p+p/p+p and 
Pb+Pb/Au+Au at ultra high energies (17.3 GeV < y^snn < 14 TeV) from the Ultra-relativistic 
Quantum Molecular Dynamics (UrQMD-v2.3) model. In addition, excitation functions of produced 
charged particle multiplicities (A^'ch) and pseudorapidity spectra are investigated up to LHC energies. 
Good agreement is observed between UrQMD and measured pseudorapidity distributions of charged 
particles up to the highest Tevatron and SppS energies. 
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High energy nucleon-nucleon and nucleus-nucleus col- 
lisions are an excellent tool to study nuclear matter un- 
der extreme conditions of temperature and density. First 
day observable is the abundance of charged particles in 
elementary (anti-)proton-proton collisions and in heavy- 
ion collisions. This allows for a first exploration of par- 
ton densities in the early stage and provides stringent 
limits for nearly all available theoretical models. It di- 
rectly reflects how much of the initial beam energy can 
be converted to new particles and it is therefore di- 
rectly linked to the stopping mechanism of the initial 
protons and nucleons. Thus, the particle multiplicity 
contains information about the entropy of the system 
and the gluon density in the first stage of the collision. 
In nucleus-nucleus collisions more particles are produced 
compared to nucleon-nucleon collisions. By scaling the 
produced particle multiplicity in Pb+Pb/Au+Au colli- 
sions by iVpart (the number of participating nucleons) it 
can be tested whether nucleus-nucleus collisions are just 
a sum of nucleon-nucleon collisions or if a more collective 
type of physics is taking place. The RMS-width of the 
charged particle pseudorapidity distribution gives infor- 
mation about the longitudinal expansion of the system. 
Starting from a model benchmark in comparison to data 
from SPS, RHIC and Tevatron, we proceed to a predic- 
tion for the charged particle density expected at LHC 
energies. 

The present letter uses the microscopic transport 
model UrQMD (Ultra-relativistic Quantum Molecular 
Dynamics) in version 2.3. This new version includes (be- 
sides other changes) a coupling to PYTHIA to allow for 
the treatment of hard pQCD interactions P, [3] . Further 
detailed explanations about the changes can be found 
in 0. Therefore just a brief introduction to UrQMD will 
be given in this paper. UrQMD is a microscopic many 
body approach to p+p, p+A and A+A interactions at 
relativistic energies and is based on the covariant propa- 
gation of color strings, constituent quarks and diquarks 
accompanied by mesonic and baryonic degrees of free- 
dom. Furthermore it includes rescattering of particles. 



the excitation and fragmentation of color strings and the 
formation and decay of hadronic resonances. Moving to 
higher energies more sub-hadronic degrees of freedom are 
available and the treatment of these is of prime impor- 
tance. In the current version of UrQMD this is taken 
into account via the introduction of a formation time for 
hadrons produced in the fragmentation of strings [3, |^, 
and by hard (pQCD) scattering via an embedding of the 
PYTHIA model. The leading hadrons of the fragment- 
ing strings contain the valence-quarks of the original ex- 
cited hadron. In UrQMD they are allowed to interact 
even during their formation time, with a reduced cross 
section where the reduction factor is defined by the ad- 
ditive quark model, thus accounting for the original va- 
lence quarks contained in that hadron P, Q . Those lead- 
ing hadrons therefore represent a simplified picture of 
the leading (di)quarks of the fragmenting string. Newly 
produced (di)quarks do, in the present model, not in- 
teract until they have coalesced into hadrons however, 
they contribute to the energy density of the system. A 
more advanced treatment of the partonic degrees of free- 
dom during the formation time ought to include soft and 
hard parton scattering 0, H, HO, 1111 and the explicit 
time-dependence of the color interaction between the ex- 
panding quantum wave-packets (l^ . However, such an 
improved treatment of the internal hadron dynamics has 
not been implemented for light quarks into the present 
model. 

The UrQMD model has been applied successfully to 
explore heavy ion reactions from BNL-AGS energies 
{Eux, = 1^-10 A GeV), over CERN-SPS energies (iJiab 
= 20A-160A GeV) up to the full BNL-RHIC energy 
(a/snn = 200 GeV). This includes detailed studies of 
thermalization ll^ . particle abundances and spec- 
tra J^, J^], strangeness production [17], photonic and 
leptonic probes J/^s [l^ and event-by-event fiuc- 
tuations j20l . [2l| . 

In the next Sections we set the stage for further inves- 
tigations by comparing UrQMD calculations with mea- 
surements performed in p+p and Pb+Pb/Au+Au colli- 
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FIG. 1: (Color online) Pseudorapidity distribution of charged 
particles in inelastic minimum bias p+p collisions from top 
SPS energies to the highest RHIC energies predicted by 
UrQMD (a). The pseudorapidity distribution of charged par- 
ticles in inelastic minimum bias p+p collisions measured by 
the UAl collaboration [S^] (b). The closed symbols indicate 
measured points, whereas the open points are reflected with 
respect to mid-pseudorapidity. The solid line represents cal- 
culations from UrQMD, in inelastic minimum bias p-|-p. 



sions from 17.3 GeV at the CERN-SPS to 1.8 TeV at 
Fermilab. This systematic comparison sets the fomida- 
tion for the following predictions for p-l-p and Pb-|-Pb 
collisions at LHC energies. 

Fig. [T] (a) shows the dNch/dr] distribution (77 being 
the pseudorapidity) for charged particles in inelastic |44| 
minimum bias p+p collisions from top SPS to top RHIC 
energies predicted from UrQMD. Fig.[l](b) presents mea- 
surements performed by the UAl collaboration 22i] for 
inelastic minimum bias p+p collisions at 540 GeV. The 
closed points show the measured region in 77, whereas the 
open points are the reflected points at 77 = 0. With in- 
creasing energy the leading hadron effect becomes more 
visible and from the gap between the humps the strength 
of the stopping effect is visible. The system is becoming 
more transparent at higher energies which is reflected 
in the change of the pseudorapidity distribution from a 



FIG. 2: (Color online) Pseudorapidity distribution of charged 
particles in inelastic minimum bias p+p collisions for differ- 
ent energies measured by the UA5 [25[ (a), CDF [ij and 
P238 (b). The closed symbols indicate measured points, 
whereas the open points are reflected with respect to mid- 
pseudorapidity. The solid line represent calculations from 
UrQMD, in inelastic minimum bias p+p. 

Gaussian to a double Gaussian shape [H, [l^l . The same 
structure is also visible for the charged particle pseudora- 
pidity distribution in inelastic minimum bias p+p colli- 
sions at y/s = 53, 200, 546 and 900 GeV measured by the 
UA5 collaboration [H (see Fig.[2](a)) and the P238 [2^ 
and CDF [27j collaboration in inelastic minimum bias 
p+p collisions at 630 GeV and 1.8 TeV collision energy 
(see Fig. [2] (b)). A difference is observed between the 
experiments P238 and CDF at 630 GeV collision energy. 
At first glance it seems that a discrepancy between the 
measurements of UAl and UA5 at 540 GeV and 546 GeV 
exists. However, in the authors assure the reader that 
both experiments agree within the error, therefore we re- 
frain from discussing possible reasons for the apparent 
discrepancies. 

The solid lines in Figs. [1] and [5] represent calcula- 
tions from UrQMD in inelastic minimum bias p+p col- 
lisions. Unfortunately, no measurements of charged par- 
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tide pseudorapidity distributions were performed for in- 
elastic minimum bias p+p collisions at SPS (17.3 GeV) 
and RHIC energies to complete the overall picture (note 
however, that pion distributions at SPS and RHIC are 
well described by the present model [3]). Comparing 
UrQMD to the measurements from the UAl (see Fig. [1] 
(b)) and UA5 (see Fig. [2] (a)) the model describes the 
UAl data on a level of w 20% and the UA5 data within 
5% accuracy. Moving to higher energies UrQMD de- 
scribes the measured peseudorapidity distribution per- 
formed by P238 (see Fig. H (b)) at 630 GeV quite well. 
Comparing UrQMD to the measurements fr'om CDF at 
630 GeV it agrees on a level of « 25%. Also here, the 
reader should notice the difference in the measurements 
between P238 and CDF at 630 GeV. For the measure- 
ments at 1.8 TeV the deviation is on the level of less than 
10 %. 

Moving on to nucleus-nucleus reactions, Fig. [3] shows 
the diV^-/dy and dNch/di] distribution in Pb-|-Pb and 
Au+Au collisions for different experiments and energies 
from SPS to RHIC energies. Fig. [3] (a) presents the 
dN/dy distribution of negativel y ch arged pions measured 
by the NA49 collaboration from 6.3 to 17.3 GeV 

(7% most central collisions for 6.3 - 12.3 GeV, 5% most 
central collisions for 17.3 GeV) center-of-mass energy. It 
is visible that UrQMD overpredicts the measurements at 
mid-rapidity by « 5% except for the ones at 17.3 GeV col- 
lision energy. Going to the higher RHIC energies (Fig. [3] 
(b)) we compare to the measurements from the PHOBOS 
collaboration 2X, 2^ ■ It is visible that the multiplic- 
ity increases with collision energy from 19.6 to 200 GeV 
(6% most central collisions). Furthermore the shape of 
the spectra is also changing as already seen for p-l-p colli- 
sions due to the fact that the the colliding nuclei become 
increasingly transparent [H, [l^l • This is reflected in the 
UrQMD prediction where the shape of the spectra is also 
changing with energy. UrQMD slightly (20%) overpre- 
dicts the measurements around mid-pseudorapidity at 
62.4 GeV and 130 GeV. 

A crucial point for particle production in A-l-A reac- 
tions is how much of the initial longitudinal motion is 
transformed to particles and transverse expansion. This 
is best characterized by an investigation of the energy (ra- 
pidity) loss of the initial nucleons. New measurements 
at SPS energies {20A - 80A GeV) [H combined with 
previously published results from AGS to RHIC ener- 
gies [13, Ha HI, [131 are available to test the predictions 
performed by the UrQMD model. Fig. 2] depicts the en- 
ergy evolution of the relative rapidity loss of the incoming 
nucleons in Au+Au/Pb+Pb reactions up to LHC ener- 
gies. The net-baryon distribution {dNg_g/dy) is made 
by using the calculated rapidity spectra for p, p, n, n, 
A, E^, S° and n~ and their anti-particles re- 

spectively. From the net-baryon distribution an average 
rapidity shift (Sy) can be calculated as follows: 
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FIG. 3: (Color online) Rapidity distribution of in Pb+Pb 
collisions at SPS energies from 6.3 to 17.3 GeV (7% most 
central collisions for 6.3 - 12.3 GeV, 5% most central collisions 
for 17.3 GeV) measured by the NA49 collaboration H, ^ 
(a). The pseudorapidity distribution of charged particles in 
Au+Au collisions (6% most central collisions, b < 3.55 fm for 
the data from 19.6 - 130 GeV and 6 < 3.65 fm for the 200 GeV 
dataset) at RHIC energies from 19.6 to 200 GeV performed 
by the PHOBOS collaboration [3^, '31', ^33, 13^ (b). The solid 
line represent calculations from UrQMD (6 < 3.9 fm for 7% 
most central Pb+Pb collisions from 6.3 - 12.3 GeV, b < 3.4 
fm for 5% most central Pb+Pb collisions at 17.3 GeV, b < 
3.6 fm for 6% most central Au+Au collisions from 19.6 - 130 
GeV and 6 < 3.7 fm for 6% most central Au+Au collisions at 
200 GeV). 

where yp is the projectile rapidity and (iVpart) the num- 
ber of participating nucleons. It is clearly visible in the 
data that {6y) /yp decreases from « 0.6 at AGS energies 
to 0.4 at top RHIC energies which indicates that the rela- 
tive baryon stopping is slightly weaker at RHIC energies 
as compared to lower AGS and SPS energies. The same 
trend is also observed in UrQMD (black line in Fig. 3]) 
where the absolute stopping follows the trend going from 
AGS to LHC energies. Another approach is also shown in 
Fig. m from a color glass condensate model [ll] (dotted 
line). In this model the authors are using the rapidity 
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FIG. 4; (Color online) Relative rapidity shift {5y}/yp as a 
function of the center-of-mass energy in relativistic heavy ion 
collisions from AGS to RHIC energies ^31,;34;, 33;, 36, 37]. The 
black line represents the prediction made by UrQMD from low 
AGS to LHC energies. The dotted line represents calculations 
from a color glass condensate model [38l ]. 



distribution of net protons {p — p) in central heavy-ion 
collisions as a testing ground for saturation physics and 
that the valance quark parton distribution is well known 
at large x, which corresponds to the forward and back- 
ward rapidity region. 

From these studies of the energy deposition (stopping) 
and particle production, we conclude that UrQMD has a 
valid basis for further extrapolations in energy and allows 
us to make predictions for LHC energies. 

The predictions for the charged particle pseudorapidity 
distributions at LHC energies are shown in Fig [5] (a) for 
inelastic minimum bias p-fp collisions at 5.5, 10 and 14 
TeV and for the 5% most central ((iVpart) 383) Pb-^Pb 
collisions at 5.5 TeV (b) (solid Hne). 

There are two complementary production mechanisms 
at LHC energies: hard parton-parton scattering and soft 
processes. Particles produced in hard scatterings are usu- 
ally created in primary collisions and are centered in a 
narrow region around mid-pseudorapidity (seen in dot- 
ted line in Fig. [5] (b)), whereas soft produced particles 
are distributed over the full pseudorapidity range (see 
dashed line in Fig. [5] (b)). At LHC energies both mech- 
anisms play an important role so that the pseudorapidiy 
distribution of charged particles shown in Fig[n](b) (solid 
line) is the sum of both processes. 

Fig. [S] (a) shows the measured number of charged par- 
ticles at mid-pseudorapidity (j^|^/j^=o) as a function of 
V^NN for p+Mcircles}j2|j2a]2^ Pb-hPb/Au-hAu 
(squares) [li,!!^ IIiniF''^'^ 



iaEd] colHsions It is 
clearly visible that in A-l-A collisions A'ch scales linearly 
with the center-of-mass energy. The difference in scaling 
with A'^part between p-l-p/p-l-p and Pb-|-Pb/Au-|-Au colli- 
sions increases with increasing center-of-mass energy. A 
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FIG. 5: (Color online) Prediction of the charged multiplicity 
pseudorapidity distribution for inelastic minimum bias p+p 
collision from ^snn = 5.5 to 14 TeV (a) and Pb+Pb colli- 
sions (5% most central collisions, b < 3.35 fm) at 5.5 TeV 
(b) collision energy from UrQMD, with PYTHIA (solid line), 
without pQCD contributions (PYTHIA) (dashed line) and for 
hard produced particles (dotted line) {b < 3.4 fm for 5% most 
central Pb+Pb collisions at 5.5 TeV). 



simple approach to extrapolate the number of charged 
particles in Pb+Pb collisions was suggested in [4l| by 
using a fit function (^|^/y=o = 0.5+0.39-ln(s)). It 
is visible that the fit function and UrQMD agree until 
top RHIC energies. At higher energies UrQMD predicts 
a higher multiplicity in central Pb+Pb collisions, espe- 
cially for top LHC energies as compared to the simple 
extrapolation. The reason for the increasing numbers of 
the multiplicity is the increase of hard collisions at LHC 
energies. When not taking hard collisions into account 
(see Fig. [S] (a)) by switching off PYTHIA and just allow 
UrQMD to have soft particle production, UrQMD would 
follow the simple linear fit function. If the LHC data fall 
on the dotted line, hard collisions are either absent at 
LHC or saturation effects do effectively suppress a large 
part of the particle production. UrQMD not only de- 
scribes the multiplicity and trend in p+p/p+p collisions 
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FIG. 6: (Color online) The energy dependence of the num- 
ber of charged particles ( '"J^'^ ) at mid-pseudorapidity for p+p 
(circles) and Pb+Pb/Au+Au (squares) collisions divided by 
.'Vpart (a). RMS width of the pseudorapidity rapidity distribu- 
tions as a function of the center-of-mass energy (b). The black 
solid line represents calculations from UrQMD for p-(-p/p4-p 
and the red solid line for Pb+Pb/ Au+ An collisions respec- 
tively. 
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FIG. 7: (Color online) Comparison of the predicted pseudo- 
rapidity distribution of charged particles (a) and the charged 
particle multiplicity at mid-pseudorapidity (b) from UrQMD 
and predictions from various other models pl| . 



(dashed line) but also in Pb+Pb/ Au+ An collision (solid 
line). Furthermore in UrQMD, if going to LHC ener- 
gies, the difference between p+p and Pb+Pb collisions 
becomes larger. 



The RMS-width [47[ is calculated by fitting the mea- 
sured pseudorapidity distribution of charged particles 
from UAl, UA5, P238 and CDF experiments for p+p 
NA50 and PHOBOS for Pb+Pb/AuAu collisions by a 
double Gaussian (see Fig. [6] (b)). An increase of 
the RMS-width is observed for p+p and Pb+Pb/ Au+Au 
collisions with the center-of-mass energy. The depen- 
dence is linear for p+p and Pb+Pb/Au+Au collisions. 
In the data, no difference between the RMS-width in 
p+p and Pb+Pb/ Au+Au is visible. UrQMD shows a 
slight difference between the RMS-width for p+p and 
Pb+Pb/ Au+Au collisions. 



To have an overall picture how the presented predic- 
tion of UrQMD compares to other approaches Fig. [7] de- 
picts the compiled results from other model predictions. 
Fig. [7] (a) shows the predicted pseudorapidity distribu- 
tions of charged particles from various models [4l| in 
comparison to UrQMD. It is visible that all transport 
models (hadronic/partonic), including UrQMD, can be 
put together in one group by predicting a similar shape 
and multiplicity. The second group are saturation mod- 
els which in general predict a lower multiplicity (also seen 
in [42]). This is also visible in Fig.[7](b) where the energy 
dependence of predicted charged particle multiplicity at 
mid-pseudorapidity is shown. At first glance it seems 
that the data would follow more the trend of a straight 
line but the major part of the models including UrQMD 
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do not favour this trend (also seen in [43|). 

In this paper we presented LHC predictions from the 
Ultra-relativistic Quantum Molecular Dynamics model 
(UrQMD). We started by showing that UrQMD describes 
the charged particle pseudorapidity spectra in p+p as 
well as for Pb+Pb/Au+Au collisions up to Tevatron 
energies. Furthermore it also describes the energy de- 
pendence of charged particles in mid-pseudorapidity in 
p+p and Pb+Pb/Au+Au collisions. The observed sim- 
ilar RMS-width of the charged particle pseudorapidity 
distribution in p+p and Pb+Pb/Au+Au collisions can 
also be described by our model. At LHC we predict 
dNch/drj |i4A ToV p-i-p « 6.3 and dNch/d-q I5.5 tcV Pb+Pb 
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